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Abstract 

The theory is developed describing transitions induced by electron impact 
between all the ground state terms of atoms and ions with configurations (2p)^ 
and (3p)^. Pull account is taken of exchange and the computed cross sections 
satisfy exactly the required unitarity bounds. A description is then given 
of the method used for the solution of the resultant coupled integro differential 
equations* 
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1 . ^ Introduction 

The interpretation of astronomical observations is helped considerably by 

a knowledge of the atomic processes which can take place in the constituent 

elements of astronomical objects. In particular, by studying the observed 

spectrum one gets some idea of their chemical composition, density and temperature. 

Many spectral lines observed in gaseous nebulae, in auroras, and in the night 

airglow are due to transitions,, among the terms of atomic systems, which involve 

no change in the electron configuration of the system. Due to the conservation 

of parity, these transitions are forbidden for electric dipole radiation but are 

allowed for electric quadrupole and magnetic dipole radiations. The transition 

probabilities are of the order (1 to 10"^)sec" , compared with 10 sec” for 

permitted dipole lines. The strongest lines in the spectra of most nebulae are 

due to the a2i ^‘ ^2*^1 ma S ne ‘ t ^ c dipole transitions in 0^^ , which has the 

configuration (Is) (2s) (2p) in both the initial and final states' . Indeed 

there is considerable astrophysical interest in all the forbidden lines arising 

from transitions between the three spectral terms of configurations" and 

q=2,3>4 i n a large number of atoms and ions. These lines can be excited by 

* 

electron impact. Consequently, it is of interest to predict theoretically the 
excitation cross sections for electrons incident on atomic systems with configurations 
(Is) 2 (2s) 2 (2p) q and (is) 2 (2s) 2 (2p) 6 (3s) 2 (3p) q . 

There is additional astrophysical interest in these transitions for at least 
two reasons. Firstly, they play an important role in determining the electron 

density distribution and the electron temperatures in the E and & regions of the 

( 2 ) 

upper atmosphere. Secondly, it has been pointed out by Branscomb and i^age}' ' 
that the absoiption of photons by atomic and molecular negative ions, especially 
C and possibly 0 , may be of considerable importance in understanding the physical 
properties of stellar photospheres. In order to compute the bound-free and free-free 
absorption coefficients it is necessary to know the radial function describing the 
motion of the electron, freed by photo-detachment, relative to the parent atom. 

This is precisely the scattering radial function. # 
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fiarly calculations of the electron induced transitions in these atoms were 
shown by Bates et al. , to exceed the conservation condition by a large factor. 
Seaton in a series of papers developed the theory of the continuum Hartree-Fock 
equations and applied it to a variety of atoms and ions. However^ owing to the 
complexity both of the algebra and of the numerical evaluation of the resulting 
equations only a few calculations of limited accuracy have been carried out. 

Recently* the theory has been extended by Seaton and applied by Shemming^^ to 
calculations of transitions in Both the methods of distorted waves and 

exact resonance, were used in the solution of the resultant coupled integro- 
differential equations. Recently, also, some exact solutions of a single second 
order integro-differential equation describing the elastic scattering on the 
ground state of C have been reported by Myerscough and UcDowell^^ and on 0 + 

by Dalgamo , Henry and Stewart N ; 

In this paper the theory of scattering of electrons by atomic systems with 
configurations 2p^ and 3q^ , where q=0 to 6, is developed. Particular emphasis 
is given to developing the formalism in a way which takes full advantage of the 
power of high speed digital computers to aid in the algebraic analysis as well 
in solving the resultant equations. It is found that if a proper treatment of 
exchange is adopted, then the post-prior discrepancy which troubled early workers 
in this field does not arise. Further, if only the terms in the ground state 
configuration of the atom or ion are coupled then the resultant integro-differential 
equation can be easily solved without further approximation using fast computers. 
There is very little point in attempting to simplify the equations to use approximate 
solutions thus introducing further unnecessary errors and ambiguities. 

It is convenient to discuss here in general terms the expected accuracy and 
thus justification of our approach for -these transitions. It is by now well known 
that the close coupling expansion can give results of dubious accuracy for excitation. 
Thusjin calculating the 1s-2s or 1s-2p excitation by electronimpact in atomic 
hydrogen the 1s-2s-2p close coupling approximation gives results much larger than 
experiment (Burke and Smith^). It is not sufficient to include just a few further 
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excited states in the expansion but the whole series must be rearranged to obtain 

faster convergence. However, for elastic scattering oh the ground Is state of 

hydrogen the approximation obtained by retaining only the Is state in the • 

expansion gives, when exchange is included correctly, ten per cent accuracy. 

This is basically because the Is state is well separated in energy from neighbouring 

. q 

states and is thus weakly coupled. In atoms and ions with configurations (2p) 
and (3p) ^ the ground state terms are also well isolated and weakly coupled to 
other levels. The transitions amongst these terns, is thus more analogous to 
elastic scattering in hydrogen than to excitation. Thus, provided all the direct 
and exchange interactions between the terns are included correctly, as in this 
paper, then the result can be expected to be of fairly high accuracy. 

In Sec. (il), the form of the trial function to be substituted into a 
variational principle is discussed. The use of the principle itself and the 
derivation of the continuum (i.e. scattering) Hartree-Fock equations is presented 
in Secs, (ill) and (IV). In the final section, (V) we report in detail the 
numerical methods we have developed to solve the coupled systems of integro- 
differential equations. 



2. The Trial Wave Function 


It ’will be assumed throughout this paper that the Hamiltonian is spin- 

S and L respectively 

independent; consequently, both total spin and orbital quantum numbers /are 
conserved* In order to take advantage of this fact we shall work in a 
representation which is diagonal in both L and S. It is well known that only 

L 

approximate solutions can be given to the problems of collision theory; Seaton 
has shown that the only consistent means of obtaining anti-symmetric wave functions 
in approximate solutions is to make the expansion explicitly antisymmetric. 

Consider the collision of an electron with an atomic system with nuclear 
charge Z and having N electrons, let the totally antisymmetric (N+l) electron 
wave function be 

\Ty ^£2*1 — • E * ( r y X ) 


where denotes the spatial and spin coordinates of the ith electron and T . =y . 

k^. 1^. S,. LSM^ Mg is the complete set of quantum numbers required to specify 

the atomic system in the state j* The atomic term is labelled L.S., the wave 

J D 

number and the orbital angular momentum of the projectile are k . and 1 . , respectively. 

0 3 

This (N+1 ) electron wave function can be expanded in terms of basis functions which 
are completely antisymmetric under interchange of the coordinates of a pair of 
target electrons 

N+1 

'r.i > -J (-I) 11 * 1 - 1 <1M )-i 2 |r ± : xl r; a> F.^r.) r’ 1 , 


...( 1 ) 


where X 1 denotes all the coordinates of the (N+1 ) electrons except those of the 

ith. and F. .(r) is the function which describes the radial motion of the impinging 

electron in the channel f. when the system was initially in the state T .* The 

^ J 

summation over in (l) is restricted in practic e to go over the terms belonging 
to the ground state conf iguration of the target. It will be shown below that the 
functions F^^r; satisfy coupled systems of second order ordinary integro-differential 
equations with the boundary conditions, 
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p iJ < r ) 


I 1 !* 1 


r — >0 


P ij (r ) V ( 6 ij sin 6 i + &ij 008 6 l) » k i >° 


~ «- lk J r ' K I i°e 2k i * 4 <0 


1 ...( 2 ) 


where = k ± r - l ± %/2 - r\ ± log 2k ± r + ...(3) 

i 

and ri i = - (z - N) / k ± . 
and oi = arg r(l i + 1 + in.) 

The radial functions F. .(r) are continuum Hartree-Fock orbitals; from the 

1 J 

properties of the surface harmonics F. . (r)will be automatically orthogonal to 

atomic discrete orbitals with orbital quantum number * 1^. For closed 

subshells, it follows from the antisymmetry cf Eq. (l) with respect to the 

interchange of any two of the electrons that there is no approximation in 

choosing F^ orthogonal to even though 1^ might equal l i# This may be 

modified slightly if we choose approximate Hartree-Fock orbitals for the t arget 

rather than exact H-F orbitals; however, in any case the overlap can be expected to 

be very small. For incomplete subshells, with 1. = 1. we can expect 

K i 

the overlap integral (F . . ,P) * 0. To take this effect into account , since it is 

1 J 

equivalent to the virtual capture of the impinging electron, we have imposed the 
condition that F. . are orthogonal to all discrete orbitals and added to Eq. (l) 
an arbitrary amount of a wave function, ¥ q , corresponding to a configuration with 
one additional electron in the incomplete subshell. Thus, the trial wave function 
for the (N+1) electron system in the state I\ initially, will be 


% ) - \ ®" (ls) 2 (2s) Z (2p) <1 : t X > + 0 i 'S Q [(ls) Z (2s) Z (2p) qVI :LS:X ] 


...( 4 ) 
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Equation (l*.) can be interpreted as a statement of configuration interaction where 
the first term represents an atomic system with one electron in a continuum 
orbit al,wtil£in the second term all electrons are in discrete orbitals. The coefficients 
C** are determined from the variational principle given in the next section . 

It is assumed that the orbitals P^(r) do hot depend upon the term (e.g. ^P, D, 

S for carbon) but only on the configuration and Z* This assumption is vital 
to our derivation and will be seen to greatly simplify the analysis* The error 
incurred by this approximation can be expected to be small • 


3# Variational Principle 


For the asymptotic normalization chosen in Eq#(2) f it can be 


( 8 ) 

shown, see Burke and Smith' ' , that 


6(L -<2/2) = 0 


...( 5 ) 


provides the Basis for a variational principle since the quantity ( L- j{/2 ) is 
stationary with respect to variations 


6F.. (r) "* 6 ($2^ coa 


...( 6 ) 


and arbitrary where 


L ji * s > *t (r i f * ) 3S 


with the Hamiltonian of the full system being 


ilT I ilT | 

Hjj I = 2 (-£ V? - 7“) + 2 r”“ 

w+ i=1 i i> j=1 I — i 


- £ jl' 


We will adopt atomic units with If = e = m = 1 throughout# 

Upon replacing the first tern of Eq#(4) by the expression given in (l) 
and substituting the result for the first in Eq# (?) we obtain 


Lja « / <S [ Z (N+l)* \ ( r iJ x" (Nf1) » % +1 ) + C k f o ]* 

A i 


*0^, -E) *,(*!> I>. 
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having used the fact that is symmetric and is antisymmetric under the 

interchange of all the coordinates of any pair of electrons. We have introduced 
the notation 


\ ( r i » 2 ^ % +1 ) = I X \ f N+1 <r N+1 > ( r N +i ) ...(10) 


N+1 


Now substitute the full form for in Eq. (9) and note that the term i=N+1 
is distinct from the other N terms to give 

1^ * f a* l 2 \ (r ± ; x" (N+1 ^ % +1 ) + c k (N+ir** o 

J % 

X [ 2 l» 1 (r J s X“ (N+1) ; 2S N+1 ) - H (r d ; x-N j^) J + (Nfl)* C 1 W o ] ...(11) 
r i 

In order to carry out the reduction of the six terms appearing in Eq # (ll) 
we consider separately, the direct term 

L iii - < *k 0,%*,). <W> *i (J- W >- -02) 

the exchange term 

L i^, » - -■< * k a. %,i). (H;E) \ (3; %> > -03) 

terms linear in the arbitrary constants C 

hk * hi ■ ^ to 1 < V <%?) $ V 1 - %„) > 

* ° k < V <%?> I *i O^,) >Q -Ov) 
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and finally terms quadratic in C 



a C k C 1 < W , (H-E) f >. 

0 O 


...( 15 ) 


To evaluate the matrix element given in Eq. (12), the functions and , 
are replaced by (10) with the composite function of target and spin-angle 
parts of the projectile wave function being expanded out using Clebsch-Gordan 
coefficients and Hartree-Fock wave functions for the target which are assumed to 
satisfy 

I ... 1 d = 2 1*“ *HF ^1 ** L i S i^ Vij^ 1 '^ ¥ HF ^1’*% ; L j S j) = 0 

...( 16 ) 

In other words, are so chosen that H^ is diagonal in that representation, 
where the N electron Hamiltonian is defined in terms of the N+1 Hamiltonian by 


1^1 = Hjj + ^ C% +1 ) + r N+1 ^ a 

where ^(x^) = - i V K ^ - 2 

r N+1 


...(17) 


Using (10), ( 16 ) and (17) in (12) we obtain 


‘J 


^N+l P ik ^ r H+1 ^ 


tfS - * 22 + k n 

' ^N+l r N+1 2 r H+1 > 




+ Y ij ^ r N+1 ^ ] P jl ^ r N+1^ ’ 


...(18) 
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where k* = 2 (E-E n (L ± sj) 


is in Rydbergs if E and Ej^ are in units of 27.2097 eV, and where 


V ij( r N+l) = & ( r i : 2 ^ K -N+1°N+1 ^ r N*1 fi> ^ r j :X r N+1 °N+1 ^ 


%•**% ^N+l ^Nfl • 

a 

This expression for V. (r) has been shown by Be2y, Tully and Regemorter to 

J 

reduce to 


'« (r) ■ 


subshells 


+ 5 S S 3q C(21 i+ 1)(2L + I)(21 + 1)(2L + 1)]^ Z (2X+1)” 1 
id J 0 X 


x (1.1, o« | Xo) (1 1 00 | Xo) W(l.L. 1. L.; L X) Z (-l) L * L i +L j +L 2 

° 1100 LgSg 

*(l L i s i/i L 2 S 2> s S <S L 2 S 2> *< ll i l ¥-2 ^ < P npV r) ’ — <’?> 


where n is the principal quantum number of the outermost incomplete p- subshell, 
(ab oojco) is a Clebsch-Gordan coefficient, W(abcd; ef) a Racah coefficient, both 
in the phase convention of Blatt, Biedenhacnand Rose^ and (q | ] S^) 

is the coefficient of fractional parentage. 

In the evaluation of the exchange term, Eq. (13), we note that the matrix 
element <H^> will include and overlap integral 

^ & h+i F ik^N+1^ ®2p^+1^ ^ •••(20) 
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■which, as discussed in Sec. (il), we are going to set equal to zero. The matrix 
element of the single electron Hamiltonian, ) , will contain a similar 

factor to (20). Consequently, ( 13 ) simply becomes 


L ik,jl = “ N 'tj+I V 1 * %+1^ - %)> ...(21) 

since those terms in Z^ with a. * N will all have overlap integrals like those 

^ 1 1 ) 

of (20). Equation (21 ) was evaluated using the methods of Hartree' ' and 

( 12 ) 

Ifocah ' and we obtain 


L Im1 ‘ l*TM P ik *ii r il 


Y l V P »V P iic P jl ‘ 21 i * ’J " 1 


n 1 

= closed 
subshells 


x ( l' Xoo jl-o) 2 - 1q [(21 i+ l)(2L i+ l)(2S i+ l)(21 j+ l)(2L 5+ l)(2S>l)F 


x Z (qL. S i U L 2 S 2 ) (qL.. JjLg Sg) i £ S ± i SS 2 ) 

L 2 S 2 

f L 2 1 L.> 

x 2 (l-l d oo| Xo) (a.1)- 1 ( 1 , 1 . oK 0 ) f , x i, K (P np p i ik P jl V» 


L i L ' 


...( 22 ) 


( a b e ; 

c d e'^is the Wigner 9 j coefficient. 

ff’ g 


' In the evaluation of the terms linear in C, both the matrix elements of 
and E contain overlap factors as in (20) and the second term in (14) becomes 

* (N+1)* 0 k Z ; dX W 0 (X) [H.j (x^) + N_ ] *i(i» % +1 ) 

1 r N+1 ,N 

- ° lc ^ T i p ii ar > 


where 


V ± (r) = (N+1 ) 2 [6^1 (q+1 LS |] L. S.) K4 ~ 2 + ^2 " f ) P np (r) 


+ Z , 2(21+1 ) V r) +3q 

n 1 

= closed 
subshells 


z 

l's' 


6 s *g(q+1 LS J] L S') 


[3(21.+ 1 ) (2L* +1 ) (2L ± +1 ) ] 


JL 

2 z 

xl 2 s 2 


(-1 ) l 2 + l +l 


+L, 


i (q L* s' |] L 2 S 2 ) 


X (q L i S i ]] L 2 S 2 ) (2X+1)" 1 (liooj 1 0) (11 ± 00 /x 0) W(1L* 1 ± L ± ; LX) 

x W(1L' 1L.; L, X) y x (? np P np ;r) P np (r) ...(24) 

It is emphasized that the terms involving C will only appear if electrons are 
scattered from target systems with incomplete outer p-sub shells. 

Finally we come to the evaluation of Eq. (15) which does not contain the 
radial functions F J Consequently, the matrix element will be simply a number 
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( 25 ) 


^ s ° k C 1 Oferf - E > 

= c k c 1 [ 2 f (q+1 LS | i l's') 2 
L S 

•'*% <r) <-* ^ * ? ■ * j <r> 
where the last matrix element is 


e^l's') -e 

+ N < ¥ I r” 1 . M If > ] 
o * N+1 ,N I o J 


n V 2(aVl) E o <Vl' 

= closed 
subshell su 


P 

np 


t * i * 

n 1 


P 

np 


) + q 6 


S S" 


2, (q+1 LS JJ L's') 

LSL" 


x (q+1 LS f] L"S' ) Z (q L* s' |j £ *>) (qL" s' |] * £) [(2^+1 )(2L"+1 )]* 

X £ .ft 


X (llool Xo) 2 (2X+1)" 1 W(1L* 1 L"; LX) W(l l' 1L" j j£ X) 


* \ < P np * <-’> 


L + L" + L + l 


...( 26 ) 


Defining D_ by 


D. . = 
id 


i/d 2 1.(1. 

= - 2 ~ 2 " -i-J 

Vdr 2 

> r 


(1.+1) 2Z ,2 




V. . + w. 0 


...(27) 


where I, . is an integral operator given by (22), 
Eq. (5) can be written as 


6 [ I /dr By ♦ 0“ S ; dr V F 

+ » J J 
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+ c 1 Z / dr V. F jk + c k c 1 (E k+1 - B) - £ ] » o, *..(28a) 

Since we intend to impose the orthogonality of the continuum function to all 
the discrete orbitals, then extra terms 

/dr P (r) F. . (r) 6-..1 j = k, 1 for all i ...(28b) 

*1 4 np ' ' li 1 » * 

where are the Lagrange undetermined multipliers, must be added to dq. (28a). 




4. Radial Equations 


Variation of (28) with respect to gives 


2 £>F D . F„ 6, +2 F. v D. 6P 6. +C V ffi 6, 

. mn mj jl kn . ik am mn In mmnln 

J X 

+ y 6F 6, + (i P (r) 6p 8, 4 &,+p P (r) 6F 6. .6, 

»m mn kn m np* mnJL 1 kn^ m np • mn 1 .1 In 

m m 


+ * 6 &mn ** 6 nl = 0 


Using Green's theorem, the second tern in (29) can be rewritten and the result 
is 

‘ * I «to [ 1 ^ \ * «! 1», W ] 

1 Til * 


• • • ( 30 ) 


+ 6. [ 2 D . F C k V + p 6. . P (r) ] ] = 0 

In fc . mi ik m nm 1 1 np ' ' 4 * 

1 m 


For arbitrary variations 6F, subject to Eq. (6), then 


] D « F » W * cl V \i w * 0 -00 


are the system of integro-differential equations for the radial functions F . 

m jl 

Variation of (28) with respect to C yields 


6 km ^ * «•** Y < p -ii + 2 / dr V, F 


i il la . 


i ik 


♦*km CJ 


- E) + C k 6^ (E n+1 - E) = 0, 


...(32) 
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which is satisfied by 


C 1 = - - E) _1 2 / dr V. (r) F^ (r). 


• • • ( 33) 


Substituting equation (33) into (31 ) gives 


Z 

5 


V) 




£ - 


^+1 


V Jr) Z I dr' V (r) F (r') 
m j J Jl 


+ 6 


V 


JJ p 
m np 


(r) = 0 


...(34) 

These coupled integro-diff erential equations are solved by the methods described 
in the next section to yield the radial functions F ^ (r). The Lagrangian 
multiplier’ s are adjusted so that the integrals (28b) are zero. The matrix 
can be simply determined from the asymptotic form of the functions F .. (r) according 
to equation (2a). From the explicit symmetry of the direct and exchange potentials 
in (34) it follows that the ^-matrix is also symmetric and thus the resultant cross 
section satisfies the required unitarity bounds. The transition matrix T is then 
defined in the usual way by the matrix relation 

T = 2 i Si/ (1- i $Q ...(35) 


and the cross section for the transition L. S. — > L . S . follows immediately 

0 0 
(1 3) 

from an analysis similar to that given by Lane and Thomas' . It is defined by 


°L.S.,L.S. = T Z Q I 2L+1) . (2S+ _ 1J 

11 i J L S * 2kT (2 L.+1)(2S.+1) 

1 , 1 , 1 1 



...(36) 


in x a Q units, where L, S, % are the total angular momentum, total spin and parity 

of the system respectively and where 1. and 1. are the orbital angular momenta of 

^ J 
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the initial and final scattered electron states respectively. 

In a typical case, for example the scattering of electrons by atomic 
oxygen, conservation of L S and x implies that (34) reduces to a set of either 
four of five coupled integro-differential equations for S=g- depending on whether 
L+x is odd or* even respectively. 


5. Computer Program to Solve Equations 

A FORTRAN program has been written to solve equations (34) for an electron 
scattered by an atom or ion of arbitrary charge and with the 2p or 3p shell 
partially filled. The program uses subroutines for the Clebsch-Gordan coefficient 
(a boo I co), the Racah coefficient W(a b c d; e f) and the Wigner 9j coefficient 


( a b e\ 

c d el and a table cf relevant fractional parentage coefficients (q L, S. I |L 5 S ? ). 
f f' g/ ^ 

With this information the program carries out the summations involved in the definitions 

of V. .(r), W. . and V. (r) and decides, for a particular L S and r. input values, 

how many equations in (34) are coupled and what the explicit form of all the 

interaction terms are. It also decides automatically on the basis of the particular 

atom or ion being considered and the energy of the incident electron, what intervals 

to use in the numerical integration of (34). The remaining information required 

by the code in order to carry out the evaluation of the, cross section defined by 

( 36 ) is the Hartree Fock orbitals for the atom or ion. It was found convenient, 

(l4) 

to read in the analytical orbitals defined by Roothaan and coworkers ' since 

they are now available for most of the atoms and ions of interest. 

The numerical solution of the resultant coupled equations (34) will now be 

described. To enable the solution of (34) to be obtained when some or all cf the 

channels are virtual it was necessary to adopt a method of inward and outward 

integration with subsequent matching to obtain a final continuous solution. 

(15) 

This is the method used by Smith and Burke v ^ . The integral terms in (34) are 

(16) 

treated non-iteratively following Harriott' ' in order that any narrow resonant 

effects may be determined without convergence difficulties (see Burke and McVicax 

(18) 

Finally the asymptotic expansion nethod of Burke and Schey ' ' was used to 

determined the ^ -matrix from the functions F^« Equation (34) can now be written 


ai*i 

dr 2 


M NE 

. \ + k l 1 “k \ (P k P k' r) P k ( r) 


+ 2GV ± 


.+ 




P 

np 


(r) 


-.(37 a) 
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where the suffix denoting the incident channel has been dropped since we require 
all possible solutions of (37) consistent with the boundary conditions (2a). 
Further the quantity C in (37) is just that given in (33 )• 

According to Hartree' 1 , the functions y^(P F; r) satisfy the following 
second order ordinary differential equation 

d^ (ryj = X (X+l) (rl) - (2X+1) P(r) F(r) 

2 * o A 

dr r r* ...(37b) 


with the boundary conditions 


r y-v (r) - r 


X +1 


r — > o 


r y x (r) ~ r 

r > « 


-X 


...(2b) 


We solve the equation (37) subject to the boundary conditions (2) 

We now introduce the following definitions: The letter M denotes the number of 
different F f s, i.e. the number of channels; NE is the total number of exchange 
terras appearing in all the F-equations* Let NA be the number of open chanrels 
(i*e. those channels for which k^ > 0). Y/e introduce the parameter NV which 
is zero if all.V^ = 0, otherwise it is unity* Let NMJ be the number of Lagrange 
multipliers in a given set of equations (37)* 

At r'=0 we define (MfNE) linearly independent solutions of the homogeneous 
system of equations i.e* 


Z" = B. Z ...(38) 

^here Z = 

is a column vector which has (M+NB) elements. Since Z(0) = 0, then (M+NE) of 
the arbitrary oonstants are fixed and the independent solutions are generated by 



- 19 - 


setting the coefficients o f the poweiv or at the origin unity one at a 

time with the others set equal to zero. These solutions are then integrated out 
with C =0. A further ( NV + N12J) independent solutions of the inhomogeneous 
system are generated by setting either C = 1 and all the » 0, or setting one of 
the equal to unity and the remainder, together with C, equal to zero. In this 
way we define NIN = M + NE + NV + NMJ linearly independent solutions at the origin. 
These solutions are then integrated out to some match point denoted r and the NIN 

v-f a 

solutions are labelled * These solutions are combined with (as yet unknown) 

coefficients where a=1 , ... , (M + NE) , C and ja to give the full solution’ over 

the range 0 ^ r ^ r^. We choose to use the Runge-Kutta method to integrate from 

r = 0 to a mesh point h ~ 0.01 and then the Numerov method from h to r . 

o 

Asymptotically, where r = r^, (M+-NA) linearly independent solutionsof Eq. (38) 
are generated by setting the coefficients of sin 0* , cos 0 # and e"* r each in 
turn to unity, and the remaining coefficients to zenq where 0'= kr - t] log 2kr and 
i*=k. The asymptotic expansion of Burke and Schey is used to calculate the function 
at and (r^- H) where H is another, stepping increment. These solutions are 
integrated in to some point r , where the exponential terms in the various 

A 

potentials might be expected to begin contributing. At r . a further NE linearly 

independent solution are generated by setting the coefficients of in turn to 

unity and all others to zero. A further (NV+NMU ) independent solutions of the 

inhomogeneous system are defined by setting C or equal to unity as in the 

inner region. The totality of outer solutions, NOUT = IT+NA + NE + NV + NMU are 

integrated from r. into r and beyond to some r '. These outer solutions are 
A O O 

labelled ^ 

In order that the solution be continuous over the whole domain 0 v< r < r D . it 
is possible to impose continuity of functions and derivatives at any point, e.g. 

^ at r Q , or the functions be continuous at a pair of points and say. The 
latter criterion was used in the computer program, and gives 2(M*NE) equations 


■ 20 - 



for the 2(M+NE) + NA + NV + NMU unknown parameters C and A further 

(NV + NMJ) equations are obtained by substituting the continuous solution given 
by Eq. (39) into ( 33 ) and also into the equations obtained by requiring (20) to be 
zero. 

The remaining NA equations are obtained by specifying that the coefficients 

of the sine part of the R.H.S. of Eq. (39) in the asymptotic region equals 

kT 2 6 . for j = 1 to NA. That is 
1 ij 

0 1 2 NA p NA+i 

^ “0 Si s "i ii + “2 h i + “na ^ i + ••• + “na + i bi + •••• 


^ i ^ NA+i since the others do not contribute 

“i * “Wi Si 


= sin 0 / + “NA+i cos ®i* ••»(40) 


where 0 / = k^r - log 2k^r . Prom Eq. (2a) 

• _i 

F^(r) k^ 2 [(sin 0^ cos 0^ - cos 0! sin $/) 6^ 


+ (cos cos 0^ + s i n sin 


- 2*1 



where 0^ = 1^ Tt/2-cr^ 'Combining this result with (l+0) gives the required 
NA equations 


cos 0 ± c# ± w - sin 0 i <u NA+i ‘ J = 2 6 


ij 


...( 41 ) 


where the superscript j denotes that the linear equations at the matching point 
for the parameters p must be solved NA times, with j characterizing these 

matched solutions. The elements of the real symmetric reactance matrix are then 
given by 

(Rij = k ± ? C sin 0 i w i J + 008 ^i “na + i J i ...( 42 ) 

Having obtained the matrix it is a simple matter of matrix manipulation 
to obtain the T -matrix and cross section using equations (35) aM (36). 


m 
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